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Abstract 



q 

^ We study deterministic discrete time exclusion type spatially heterogeneous par- 

■ ticle processes in continuum. A typical example of this sort is a traffic flow model 
with obstacles: traffic lights, speed bumps, spatially varying local velocities etc. 
Ergodic averages of particle velocities are obtained and their connections to other 

■ statistical quantities, in particular to particle and obstacles densities (the so called 
Fundamental Diagram) is analyzed rigorously. The main technical tool is a "dy- 
namical" coupling construction applied in a nonstandard fashion: instead of proving 

£\j ' the existence of the successful coupling (which even might not hold) we use its pres- 

ence/absence as an important diagnostic tool. 

in 
o 

1 Introduction 

In 1970 Frank Spitzer introduced the (now classical) simple exclusion process as a Markov 
chain that describes nearest-neighbor random walks of a collection of particles on the one- 
dimensional infinite 1 integer lattice. Particles interact through the hard core exclusion 
rule, which means that at most one particle is allowed at each site. This seemingly very 
particular process appears naturally in a very broad list of scientific fields starting from 
various models of traffic flows [14, 9, 7, 2, 3], molecular motors and protein synthesis in 
biology, surface growth or percolation processes in physics (see [15, 5] for a review), and 
up to the analysis of Young diagrams in Representation Theory [6]. 

From the point of view of the order of particle interactions there are two principally 
different types of exclusion processes: with synchronous and asynchronous updating rules. 
In the latter case at each moment of time a.s. at most one particle may move and 
hence only a single interaction may take place. This is the main model considered in 
the mathematical literature (see e.g. [13] for a general account and [1, 8, 10] for recent 
results), and indeed, the assumption about the asynchronous updating is quite natural 
in the continuous time setting. The synchronous updating means that all particles are 
trying to move simultaneously and hence an arbitrary large (and even infinite) number 
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of interactions may occur at the same time. This makes the analysis of the synchronous 
updating case much more difficult, but this is what happens in the discrete time case. 2 
This case is much less studied, but still there are a few results describing ergodic properties 
of such processes [2, 3, 7, 9, 14]. 

Recently in [4] we have introduced and studied the synchronous updating version of 
the simplest exclusion process (TASEP) in continuum. 3 Now our aim is to extend these 
results to the case when the media is heterogeneous, i.e. contains obstacles. The idea is 
as follows. Consider a (countable) collection of particles performing (random) walks on 
a real line with hard core type interactions (we assume that the particles cannot outran 
each other). Assume also that on the real line there is a (countable) number of obstacles. 
To overcome an obstacle a particle needs to spend some additional time on it. Thus we 
have two types of interactions: between particles and with static obstacles. The principal 
novelty here is that the presence of obstacles leads to a spontaneous creation of "traffic 
jams" (particles clusters) near obstacles. Indeed, it is easy to construct initial particle 
configurations having no clusters of particles but such that these clusters will be created in 
front of obstacles under dynamics. This feature was completely absent in the case without 
obstacles: a cluster may only disappear but never be created. We start the formalization 
with the simplest version of the model and later in Sections 6,7 consider what happens in 
a more complicated setting. 

By a configuration x := {xi\ ie i we mean an ordered (i.e. Xi < Xj+i Vi G Z) bi-infinite 
sequence of real numbers Xi G KL The union X of all such sequences plays the role of 
the phase space of the system under consideration. Consider also a special configuration 
z G X such that Zi < Zj+i Vi G Z. We refer to elements of x as to positions of particles 
and to elements of z as to positions of obstacles. 
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Figure 1: TASEP in heterogeneous continuum. 
Let A, = Ai(x) := x i+ i — X{ and let 

Aj = Ai(x, z) : = min ^x i+i - x h mm(zj : Zj > Xi) - x^j 

stands for the minimum between the distances from the point x\ to the point x\ +1 and to 
the next obstacle (see Fig. 1). We refer to Aj, Aj as gaps and modified gaps corresponding 
to the i-th particle in x. 

For a given real v > (to which we refer as a maximal local velocity) and the config- 
uration z we define a map T : X — > X as follows: 

(Tx)i := Xi + min (Ai(x,z),v) Vi G Z. (1.1) 

2 if one do not consider some "artificial" updating rules like a sequential or random updating. 
3 Some other continuous space generalizations were considered e.g. in [12, 16]. 



The dynamical system (T, X) describes the collective motion of particles discussed above. 
If Aj(x, z) < v we say that the i-th particle is blocked (meaning that its motion is blocked 
by the % + 1-th particle or an obstacle) at time t and free otherwise. By a cluster of 
particles in a configuration x G X we mean consecutive particles with gaps Aj(x) < i> 
having no obstacles between them. 

Associating the number of iterations of the map T with time we often use the notation 
T l x = x t := {x\}. In this terms the quantity £' := min ^Aj(x*, 2), f j plays the role of the 

local velocity of the i-th particle in the configuration x = x° ai time t > and thus the 
dynamics can be rewritten as 

xt +1 :=xi + & (1.2) 

A more general setting including varying/ random velocities and waiting times will be 
considered in Sections 6,7. 

It is of interest that even without obstacles (i.e. if z — 0) the behavior of the deter- 
ministic dynamical system (T, X) is far from being trivial. In [4] it was shown that this 
system is chaotic in the sense that its topological entropy is positive (and even infinite). 

Our main results are concerned with the so called Fundamental Diagram describing 
the connection between average particle velocities V(x) and particle/obstacle densities 
p(x),p(z) (see Section 2 for definitions) and technically the analysis is based on a (some- 
what unusual) "dynamical" coupling construction (see Section 3 and also [4]). 

For a given v > and a configuration of obstacles z denote by z the extended con- 
figuration of obstacles obtained by inserting between each pair of entries Zi, Zj + i new 
L(^i+i — z i)/ v \ 'virtual' obstacles at distances v between them starting from the point Z{. 
Here \u\ stands for the integer part of the number u. 

Theorem 1 (Fundamental Diagram) Let p(x),p(z) be well defined. Then 

V(x) = mm{l/p(z), l/p(x)}. (1.3) 

Therefore the phase space is divided into two parts: gaseous {(x,z) : p(x) < p(z)} (con- 
sisted of configurations of eventually non-interacting particles) and liquid {(x, z) : p(x) > 
p(z)} (where clusters of particles are present at any time). 




Figure 2: (left) Phase Diagram, (right) Fundamental Diagram. 

These results are illustrated in Fig 2. It is worth note that a naive argument tells 
that the average velocity of a particle in an infinite system is essentially the minimum 
of the average inter-particle distance and the average inter-obstacle distance. Theorem 1 



shows that this is absolutely not the case: instead of the average inter-obstacle distance 
one needs to take into account the average distance between the elements of the extended 
configuration z. As we shall see the latter is very different from the former, e.g. we always 
have p{z) > l/v independently on p(z) (which even might be not well defined). Actually 
the mere fact that the correspondence between the average velocities and particle/obstacle 
densities is one-to-one comes as a surprise especially without any regularity assumptions 
for the positions of obstacles. 

The main steps of the proof are as follows. First we show that for given configurations 
of particles and obstacles upper/lower average particle velocities are the same for all 
particles. Now to compare average velocities in different particle configurations with the 
same configuration of obstacles we develop a special dynamical coupling construction. 
Applying this construction we prove that that upper/lower average particle velocities are 
the same provided the same particle densities. Thus to calculate the dependence of average 
velocities on densities it is enough to construct for each density a single configuration 
having this density for which we are able to calculate its velocity explicitly. To this end 
we consider an auxiliary zero-range lattice process whose trajectories are invariant under 
the dynamics of our original system without obstacles. To calculate the average velocity 
in the zero-range lattice process one uses corresponding results obtained in [4]. 

Despite that various couplings are widely used in the analysis of Interacting Particle 
Systems (IPS) (see e.g. [13]), applications of our approach is very different from usual. In 
particular, we do not prove the existence of the so called successful coupling (which even 
might not hold) but instead use its presence/absence as an important diagnostic tool. 
Remark also that typically one uses the coupling argument to prove the uniqueness of the 
invariant measure and to derive later other results from this fact. In our case there might 
be a very large number of ergodic invariant measures or no invariant measures at all (e.g. 
the trivial example of a single particle performing a skewed random walk). This example 
reminds about another important statistical quantity - average particles velocity. The 
dynamical coupling will be used directly to find connections between the average particle 
velocities and other statistical features of the systems under consideration, in particular 
with the corresponding particle densities. 

It is worth note that all approaches used to study discrete time lattice versions of 
IPS are heavily based on the combinatorial structure of particle configurations. This 
structure has no counterparts in the continuum setting under consideration. In particular 
the particle - vacancy symmetry is no longer applicable in our case. This explains the 
need to develop a fundamentally new techniques for the analysis of IPS in continuum. 
The presence of obstacles also prevents the direct application of the scheme developed in 
[4] for the spatially homogeneous case. 



2 Basic properties 

Here we shall study questions related to particle densities and velocities. 
By the density 

,(*,/):= #{ieZ m *' £/} 

of a configuration x G X in a bounded segment I — [a,b] eR we mean the number of 
particles from x whose centers Xi belong to / divided by the Lebesgue measure |/| > of 



the segment I. If for any sequence of nested bounded segments {/„} with \I n \ n —¥ oo the 
limit 

p(x) : = lim p(x,I n ) 

is well defined we call it the density of the configuration x G X. Otherwise one considers 
upper and lower (with respect to all possible collections of nested intervals /„) particle 
densities p±(x). 

Remark 1 If 3p(x) < oo then \x n — x m \/\n — m\ \ j p(x). 

Lemma 2 The upper /lower densities p±(x l ) are preserved by dynamics, i.e. p±(x l ) = 
p±(x t+1 ) Vt>0. 

Proof. For a given segment I G R the number of particles from the configuration x l G X 
which can leave it during the next time step cannot exceed 1 and the number of particles 
which can enter this segment also cannot exceed 1. Thus the total change of the number 
of particles in / cannot exceed 1, because if a particle leaves the segment through one of 
its ends no other particle can enter through this end. Therefore 

|p(xM)- P (^ +1 ,/)|.|/|<i 

which implies the claim. □ 



By the (average) velocity of the i-th particle in the configuration iGlat time t > 
we mean ^ 

V(x,i,t) :=J^ = {x\-x°)/t. 

1 8 = 

If the limit 

V(x,i) := lim V(x,i,t) 

t^oo 

is well defined we call it the (average) velocity of the i-th particle. Otherwise one considers 
upper and lower particle velocities V±(x,i). 

Lemma 3 Let i£l and let p(z) < oo be well defined. Then \ V(x,j, t)—V(x, i, t)\ < -^> 
a.s. G Z. 

Proof. It is enough to prove this result for j — % + 1. Consider the difference between 
(average) velocities of consecutive particles 

V(x,i + l,t)-V(x,i,t) = ^±L_^±i_^_^ 

t ' t 



The last term vanishes as t — > oo and it is enough to show that the same happens with 
A*/i. Here and in the sequel we use the notation A- = Aj(x*). 



Normally in the deterministic setting gaps between particles A* are uniformly bounded 
(see e.g. [2, 3, 4]). Surprisingly in the present setting this is not the case and we are only 
able to show that A\ may grow not faster than o(t), which fortunately is enough for our 
aims. 

To prove this estimate we introduce a new configuration having only two particles 
x := {xi, x 2 } which have the same initial positions as the i-th and the (i + l)-th particles 
of x, i.e. X\ = x®,x 2 = Denoting x l := T l x and A* := x\ — x\ we observe that 

A* > A* Vt > 0. 

Consider the movement of the leading particle in the process x. There is the first 
moment of time t 2 '■= t 2 (x 2 ) when this particle encounters an obstacle from z. Denote by 
Zk the position of this obstacle in the extended configuration z. By the definition of the 
dynamics T l z k = z k+t Vt > 0. 

The movement of the trailing particle in the process x is a bit more complicated since 
additionally to obstacles it may be stopped by the leading particle. Anyway there is 
the first moment of time ti := ti(xi,x 2 ) > t 2 when this particle encounters the obstacle 
located at z k . For each t > t\ these two particles move synchronously along the "chain" z. 
Thus the growth of A* is completely determined by the concentration of entries in z (i.e. 
by p(z)). Due to the assumption about the existence of the density of z these fluctuations 
cannot exceed o(t) and hence A*/t < A*/t < -^> 0. □ 



Corollary 4 Under the assumptions of Lemma 3 the upper and lower particle velocities 
V±(x,i) do not depend on i. 

Remark 5 (a) If the density p(z) is not well defined, the distance A* may grow linearly 
with time and the limit average velocities might differ for different particles, (b) The 
existence of p(z) < 00 does not imply the existence of p(z) < 00 but fir "typical" z this 
implication takes place (see Section 7). 

Let us estimate upper/lower densities of the extended configuration z of obstacles for 
a given configuration z with p(z) > and a given local velocity v > 0. 

Lemma 6 max{l/i>, p(z)} < p~(z) < p+{z) < p(z) + l/v. 

Proof. Observe that for each i £ Z the spatial segment [zi, ^+1] contains exactly 
\_(zi+i — Zi)/v\ G [(z i+ i — Zi)/v — 1, (zi+i — Zi)/v] elements. Therefore Vn G Z + 

n— 1 



p(z,[z ,z n }) = in + ^2l(z i+ i - Zi)/v\ )/(z n -z ) 



i=0 

'n-l 



= p(z, [z ,z n ]) + \^2l( z i+i - Zi)/v\j /{z n ~ Zq) 
< p{z, [z ,z n ]) + 1/v. 

Similarly 

p(z, [z ,z n }) > p(z, [z ,z n ]) - n/{z n - Zq) + l/v = l/v. 

Passing to the limit as n — > 00 and using that z contains all elements of z we get the 
result. □ 



3 Dynamical coupling through particle's overtaking 



Consider two independent particle processes By i{x f ) we denote the i-th particle 

of the process x* (i.e. x\ is the location of the particle i{x l ) at time t). 

We say that the particle i(x t ) overtakes at time t > the particle if x i _1 < ^f 1 
and x\ > Xp and denote this event as i(x l ) > j(x f ). 

Lemma 7 If i(x t ) ^ j(x f ) then 
(a) flneZ: n(x') ^ 

/3m G Z : j(^) ^ m ( xt )- 
Additionally i(x l ) ^ anc ^ (j ~~ might happen only if x* = x t j_ 1 . 

Note that there is no contradiction between two parts of this Lemma since the former 
part concerns a particle being overtaken while the latter concerns a overtaking particle. 

Proof. By the construction of the process under consideration and the definition of the 
overtaking we have: 

x i-l ^ x i < x j 5: &j < X i < (3-1) 

Now (a) follows from the observation that by (3.1) neither particles preceding nor 
particles succeeding it may overtake the particle i(x l ) at time t. 

If (b) would hold then x'f 1 < x l ~ l <x l m < i*. Thus by (3.1) we get f* -1 < x l ~ l < x p 
which might be possible only if x l f~ l = x 1 -. The latter contradicts to the assumption about 
the overtaking. 

To show that the event discussed in the additional part may take place, consider a 
pair of configurations satisfying for some i,j,t G Z the following inequalities: 

x]-_\ < x\r x < x'r 1 = x]r + \ = x\~\, max (a^" 1 , z), A^x 1 ' 1 , zj) < v. 

Then *• which implies that i(x t ) °-> j(x t ) and (j — l)(x*) c ->- i(x t ). Assume 

now that on the contrary x\ ^ If x\ > x^_ x then x^_ x < < x\ which contradicts 

to (j — l)(x l ) i(x t ). Similarly x\ < x l j_ x implies x\ < x^_ x < which contradicts to 
i(x t ) <-» j(x l ). □ 



Let us introduce the dynamical coupling of the processes x t ,x t which consists in a 
consequent pairing of particles of opposite processes. The pairing in our deterministic 
setting is a pure formal action which does not change the dynamics. The idea is that if a 
particle overtakes some particles from the opposite process it becomes paired with one of 
them. The construction starts at time t — and initially all particles are assumed to be 
unpaired. 

Denote by 

■hir 1 ): {./( Z : itf) ^ 
the set of particles overtaken simultaneously by the particle i(x t ) at time t, by 

i,(x t ) := minjoo, min{j G Ji(x l ) : ^(f'" 1 ) is paired}} 

the paired particle with the minimal index among them, and finally 




- 1 if i,(x l ) - 1 G Ji{x l ) 
oo otherwise 



In words, i,(x l ) stands for the paired particle in J\ having the minimal index, and 
i {x t ) stands for the unpaired particle in Ji{x l ) having the maximal index among those 
preceding i 9 (x l ). A typical connection between those indices and positions of the i-th. 
particle at times t — 1 and t is shown in Fig 3. 



x -e- 



Figure 3: Connections between indices and particle's positions. 

To simplify the description of the dynamical coupling we shall use a diagrammatic 
representation for coupled configurations, where paired particles are denoted by black 
circles and unpaired ones by open circles, and use the upper line of the diagram for the 
x-particles (i.e. particles from the x-process) and the lower line for the x-particles. In 
this representation a typical pairing event looks as follows 



oo 

oo 



(3.2) 



Here =>- corresponds to the dynamics and — > to the pairing, while *, * stand for two 
different mutual pairs of particles created after the particle overtaking under dynamics. 

Now we are ready to define the pairing rigorously. We proceed first with all re-particles, 
overtaking some i-particles at time t, and then with all f-particles, overtaking some x- 
particles at this time. 

Let the overtaking takes place for i-th x-particle at time t, then if 

(a) i(x t ^ 1 ) is paired and i (^*) < oo we re-pair these particles: 

• *- , • * , • * 

• 00*0 • oo* o ' o o»* o 

(b) i(x l ~ l ) is unpaired and 

(b') < oo and x\ > Xj.( x t) we re-pair these particles: 

o •* , o» * »o * 

o '* o» * ' o» * 

(b") else if ? (^') < oo this unpaired particle forms a new pair with 

O ★ , O * • it 

oo* '' oo * ' o» * 

The pairing rules when the overtaking takes place for i-th x-particle are exactly the 
same except for the exchange of x by x and vice versa. 

The complexity of these rules reflects that first a particle may overtake simultaneously 
several particles from another process, and second an arbitrary number of particles may 
share the same spatial position. In particular, in the following event we have: 



rather than — > • °. Indeed, the trailing x-particle overtakes both two trailing 
x-particles and since the re-pairing rules in (b) demand strict inequalities only the pairing 
of trailing unpaired particles takes place. 

By Lemma 7 each overtaking event may be treated separately. Therefore we apply 
the pairing rules first for all x-particles, overtaking some x-particles, and then for all 
x-particles, overtaking some x-particles. 

According to the definition, particles from the same pair move synchronously until 
either the movement of one of them is blocked or one of the members of the pair is 
swapped with an unpaired particle from the same process. It is convenient to think about 
the coupled process as a "gas" of single (unpaired) particles and "dumbbells" (pairs). 
A previously paired particle may inherit the role of the unpaired one from one of its 
neighbors. In order to keep track of positions of unpaired particles we shall refer to them 
as x- and x- defects depending on the process they belong. 

We say that a pair of configurations (x, x) is proper if for each two mutually paired 
particles the open segment between them cannot exceed v, does not contain neither a 
defect nor an obstacle (i.e. the situations ,° *, ,|* cannot happen), and there are no 
crossing mutually paired pairs (**.*). 

Lemma 8 If the pair of configurations (x* -1 , x t_1 ) is proper for some t > 1, then the pair 
(x*,x*) is proper as well. 

Proof. The situation ,° * may happen if either a defect overtakes the trailing particle 
in a pair, or if a pair is born around this defect. Both such possibilities contradict to the 
definition of pairing. 

Assume that at time t there is a pair of mutually paired particles i(x t ),j(x t ) such that 
the open segment between them contains an obstacle: x\ < Zk < x*-. Assume also that 
this pair was present at time t — 1 as well. According to the construction of the dynamics 
this may happen only if x- -1 < z^ — v and Zk = x* _1 , which contradicts to the assumption 
that the pair (x* _1 ,x* _1 ) is proper: namely the distance between members of the same 
pair exceeds v. On the other hand if the particles i(x t ),j(x t ) were not paired at time t — 1 
and the pair is just created at time t then one of the particles should overtake another at 
this time. Hence the distance between these particles cannot exceed v. 

It remains to show that the last property still holds for "old" pairs. In order to enlarge 
the distance between the existing mutually paired particles one of them should be blocked 
by another particle or by an obstacle. On the other hand, the obstacle cannot belong to 
the open segment between pair members and hence it may block only the leading particle 
in the pair, which may only decrease the distance. 

The blocking particle might be paired or unpaired. The former case implies that the 
'companion' of the blocking particle is at distance at most v and hence it will block the 
enlargement of the distance by more than v. In the letter case the non-blocked paired 
particle overtakes the unpaired one and hence the pair under consideration will be re- 
paired. 

The observation that in the moment of the creation of a pair the distance between its 
members cannot exceed v finishes this part of the proof. 

The analysis of the absence of crossing mutually paired pairs is completely similar to 
the absence of defects between elements of a pair and therefore we skip this point. □ 



Denote by p u (x, I) the density of the x-defects belonging to a finite segment /, and by 
p u (x) := p u (x,M.) the upper limit of p u (x,I n ) taken over all possible collections of nested 
finite segments I n whose lengths go to infinity. 

We say that a coupling of two Markov particle processes x*,x* is nearly successful if 
the upper density of the x-defects p u (x) vanishes with time a.s. This definition differs 
significantly from the conventional definition of the successful coupling (see e.g. [13]), 
which basically means that the coupled processes converge to each other in finite time. 

Lemma 9 Let x,x G X with p(x) = p(x) > 0, and let there exist a nearly successful 
dynamical coupling (x*,x*). Then 

\V(x,0,t) -V(x,0,t)\ ^0. 



at time t we set n t 



Proof. Consider an integer valued function n t which is equal to the index of the x- 
particle paired at time t > with the 0-th re-particle. If the 0-th x-particle is not paired 

\n t -i if t > 
[0 if t = ' 

To estimate the growth rate of \nt\ at large t observe that n t changes its value only 
at those moments of time when the 0-th x-particle meets a x-defect. By the assumption 
about the nearly successful coupling at time t >> 1 the average distance between the 
defects at time t is of order 1/ p u (x l ) while the amount of time needed for two particles 
separated by the distance L to meet cannot be smaller than L/(2v). Therefore the fre- 
quency of interactions of the 0-th x-particle with x-defects may be estimated from above 
by the quantity of order p u (x*) — ^> 0, which implies n t /t — ^> 0. 

Now we are ready to prove the main claim. 



\V(x,0,t)-V(x,0,t)\ = \(xl-x° )-( 

< 14-41/* + 
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\n t \ + \x° -x° \/t. 



The 2nd addend is a product of two terms \n t \/t and |x* — 4I/I n *l- As we have shown, 
the 1st of them vanishes with time. If \n t \ is uniformly bounded, then the 2nd term is 
obviously uniformly bounded on t. Otherwise, for large \n t \ by Remark 1 and the density 
preservation the 2nd term is of order l/p(x), which proves its uniform boundedness as 
well. Thus the 2nd addend goes to as t — > oo. Noting finally that the 1st and the last 
addend also vanishes with time at rate 1/t we are getting the result. □ 



Lemma 10 Let p(x) = p(x) and let in the coupled process 3 a (random) moment 

of time tij < oo such that x\ > x* for each t > ty. Then the coupling is nearly successful. 

Proof. By the assumption each x-particle will overtake eventually each x-particle located 
originally to the right from its own position and thus will form a pair with it or with one 
of its neighbors (if they are so close that were overtaken simultaneously). Thus the 
creation of pairs is unavoidable. To show that the upper density of defects cannot remain 
positive, consider how the defects move under our assumptions. Assume that at time 



t > the i-th x-particle is paired with the j-th i-particle. Then by Lemma 8 in order to 
overtake at time s > t the j-th x-particle significantly (by a distance larger than v) the 
i-th re-particle necessarily needs to break the pairing with the j-th x-particle. Thus by 
the definition of the dynamical coupling either a x-defect overtakes the j-th x-particle: 
°. * — > ,° * — > .* °, or the i-th x-particle overtakes a i-defect: . * D — > . * — > *. 
(Otherwise this pair will not be broken.) Therefore during this process the x-defects move 
to the right while the x-defects move to the left. Hence they inevitably meet each other 
and "annihilate". The assumption about the equality of particle densities implies the 
result. □ 



4 Auxiliary lattice zero-range process 

Consider now a lattice process which we shall need in the sequel. This process is defined 
on an integer lattice Z occupied by a bi-infinite configuration of particles y ordered with 
respect to their positions y^ i.e. yi < y^\ Vi. The dynamics is defined as follows. 
For each i e Z consider all particles occupied the site i and choosing the one with the 
largest index (say fcj) among them we move this particle by one position to the right, i.e. 
Uki '■= Uki + ^- This is a deterministic version of the so called zero-range process on Z with 
parallel updating rules illustrated in Fig 4. 




i — 1 i i + l 



Figure 4: Zero-range process. 

Setting v — 1, z = and assuming that X{ G Z Vi 6 Z we observe that in this case 
(without obstacles) the trajectory T l x coincides with the trajectory of the zero-range 
process. Therefore according to [4] the average particle velocity for the zero-range process 
is equal to 

V(y) = min{l, l/p(y)}. (4.1) 

Assume now that the lattice under consideration is not uniform (e.g. Z) but the 
distance between the i-th and (i + l)-th site is equal to a nonnegative number for each 
% G Z. Considering the zero-range process on this heterogeneous lattice, but assuming that 
the site's occupation is the same in both cases we are able to calculate the corresponding 
statistical quantities. 

Denote by y a configuration on a heterogeneous lattice Z in which distances between 
the i-th site and the (i + l)-th one are given by the sequence of numbers {U}- Let y be 
a configuration on Z such that the i-th particle of the configuration y is located on the 
site whose index coincides with the index of the site occupied by the i-th particle of the 
configuration y. 



Lemma 11 The deterministic zero-range processes on Z and Z starting with configura- 
tions y, y described above for each t > preserves the connection between the particle 
configurations. Therefore the particle velocities in these processes satisfy the relation 

V(y) = V(y)/p(Z), (4.2) 

where p(Z) is assumed to be positive and is defined as the density of the particle configu- 
ration having exactly one particle at each site of Z. 

Proof. The first claim follows from the definition of the zero-range process. Denote 
L(t) := yl - y° Vt. Then 

(*$ - fo)/t = ] £ ^ = Tzk£ e M ^ v{v)lp{i) - 

3=0 ^ ' 3=0 

□ 



5 Proof of Theorem 1 

Lemma 12 Let x, x G X p{x) = p{x) and let V(x) be well defined for given v,z. Then 
V(x) is also well defined and V(x) = V(x). 

Proof. Let x, x e X p := {z G X : p(z) = p} be two admissible configurations of 
the same particle density. If one assumes that the dynamical coupling procedure leads 
to the nearly successful coupling of particles in these configurations then by Lemma 8 
the assumptions of Lemma 9 are satisfied and hence \V(x,0,t) — V(x,0, t)\ which 
by Lemma 3 implies the claim. In general the assumption about the nearly successful 
coupling may not hold, 4 however as we demonstrate below the pairing construction is still 
useful. 

Define random variables 

: =x*-4 ijeZ, tez . 

Then 

V(x,t,t)-V(x,j,t) = Wyt-W°/t. 

Since by Lemma 3 the differences between average velocities of different particles belonging 
to the same configuration vanish with time it is enough to consider only the case % — j — 0. 
For Wq there might be three possibilities which we study separately: 

(a) lim W l m /t = 0. Then \V(x,0,t) - V(x,0,t)\ < \W^\/t + \W$ \/t *-=^> 0, which by 

t— >oo 

Corollary 4 implies that the sets of limit points of the average velocities coincide. 

4 Consider e.g. the setting with 1/p > 5v and the configurations Xi := i/p and Xi := i/p + 2v, and 
assume that there are no obstacles. Then p(x) = p(x) = p, V(x) = V(x) = v but no pair will be created. 



(b) limsup Woo/t > 0. Then Vi 6 Z the i-th particle of the x-process will overtake 

eventually each particle of the i-process located at time t — to the right from the 
point x®. This together with the assumption of the equality of particle densities 
allows to apply Lemma 10 according to which the coupling is nearly successful. On 
the other hand, by Lemma 8 the distance between mutually paired particles cannot 
exceed v. Therefore by Lemma 9 we have \V(x,0,t) — V(x, 0,t)\ *— 0, which 
contradicts to the assumption (b). 

(c) limsupW^Q /t < 0. Changing the roles of the processes one reduces this case 

t— >oo 

to the case (b). 

Thus only the case (a) may take place. □ 

To apply this result to prove Theorem 1 one needs to construct for each particle 
density a a configuration having this density for which we are able to calculate its velocity 
explicitly. There are two possibilities to realize this idea. In both cases we use the auxiliary 
lattice zero-range process y constructed in Section 4 on the heterogeneous lattice Z whose 
i-th site coincides with the location of the i-th element in the extended configuration of 
obstacles z. 

The first construction is as follows. Choose an arbitrary initial configuration of parti- 
cles y of a given density for the zero-range process. By Lemma 11 and the relation (4.1) 
we obtain the relation (1.3). Observing that the trajectory of the zero-range process y l 
coincides with the trajectory of our original process T l y we get the desired result. 

An alternative way to derive the relation (1.3) is to construct a specific initial con- 
figuration of the zero-range process, for which we can calculate the corresponding av- 
erage velocity directly. The key observation here is that for a: := z obviously we have 
V{x) = l/p{~z). 

For a given x £ X for which 3 < p(x) < oo and a G 1+ we define the configuration 
ax 6 X in three steps depending on arithmetic properties of a: 

(a) kx for a = k G Z + ; 

(b) ^x for a = ^, k = mn + k' and m, k! G Z , k' < n G Z + ; 

(c) ^x — y ax if ^ — y a. 

Here by kx we mean a configuration in which each particle in x is replaced by k 
particles. To get y := ^x we divide x into blocks having n consecutive particles each and 
construct y := {yi} as follows. For each % at location xi we set exactly m particles at 
positions yj = X{. Besides at positions of each of the first k' particles in each of the blocks 
we add an additional particle and enumerate the particles according to their positions. 

If a ^ k/n one considers a sequence of its rational approximations: k n /n n — ^> a. For 
each n according to the rule (b) we construct a configuration i'""' := ^x. Observe that 
for £ > 1 the configuration x^™^ differs from x*- ™-* only by the presence of some additional 
particles in each of the blocks of length n£. Using the existence of p(x) one proves that 
the limit configuration exists, and we denote the latter by ax. 

Now for any a G K+ we construct x := az and calculate V(az) as follows. We start 
with the case a — 1, i.e. x — z. In this case obviously Tx coincides with the left shift of 
the configuration x. Therefore (z t — zq) • p(z, [So, z t ]) = t Wt > and thus 

V(x,0,t) = (z t -z )/t = l/p(z, [z ,z t ]) ^ l/p(z). 



If a — | the configuration x := a z is "almost" spatially periodic: it consists of 
blocks of the configuration z having equal number of n particles. The spatial periodicity 
immediately implies the existence of the average velocity V(x). Observe now that for 
any rational a < 1 the particles in the configuration az on average move exactly as 
in the original configuration z (i.e. V(az) = V(z)). If a — ^ > 1 then the direct 
calculation gives V(x) = ^V(z). Passing to the limit as ^ a we get eventually 
'-k, if a < 1 
otherwise 



V(a2) := < ^ 7 . Here one uses that > V(±z) W G Z^ 



ap(z) 



6 Varying waiting times and local velocities 

In this Section we consider a more general model which includes both varying waiting 
times and local velocities. Namely we assume that the waiting time at the j'-th obstacle 
is equal to Tj G {1,2, .. . ,r}, and the local velocity for particles in the spatial segment 
(zj,Zj+i) is equal to Vj G (0, v]. It is useful to think about this model as a road divided 
into parts separated by obstacles (e.g. traffic lights) and having different quality of the 
pavement, which we describe by varying from part to part local velocities. 

If Tj > 1 it might be possible that a new particle is coming to an obstacle when 
some preceding particles are still waiting there. We assume that for the given particle 
the waiting time Tj at the j-th obstacle starts only when the succeeding particle leaves 
it. The movement of a particle waiting at a certain obstacle in the next moment of time 
depends on the time which this and preceding particles already spent there. Therefore 
this model is non Markovian. To cure this pathology using methods developed in [3] (for 
a very different situation) one adds a new variable type (for each particle) whose value is 
equal to the amount of time the particle will wait at an obstacle if it is located on the 
obstacle and zero otherwise. 

Strictly speaking in order to use the notion of the coupling this is important (since it 
is defined for a pair of Markov processes). Nevertheless in the deterministic setting which 
we consider in this paper this issue is not crucial (as we shall see) and we proceed without 
this Markovian extension. 

Denote by t,v the collections of waiting times and local velocities corresponding to 
corresponding obstacles. If Tj = r and Vj = v we recover the previous setting. Surprisingly 
the analysis of this general setting is very similar to the previous one. Therefore we shall 
discuss only the points when the analysis differs. 

First we refine the configuration of obstacles z exchanging the original i-th obstacle by 
Ti + 1 consequent obstacles located at the same position Zj but having waiting times, i.e. 
leaving one these new obstacles the particle immediately moves to the next of them. The 
resulted ordered collection of obstacles we again call z. To take care about the change of 
indices we change also the collection of local velocities inserting r, unit velocities (before 
the original element Vi) corresponding to new obstacles and re-enumerating them. 

Obviously for any configuration of particles x the movements of its elements in the 
case of the original configurations of obstacles and velocities and refined ones are exactly 
the same. 

To construct the extended configuration of obstacles z we proceed as follows: between 
the obstacles Zi and Zi+i we insert [(z i+ i — Zi)/vi\ virtual obstacles at distance Vi starting 
from the point zi for each i e Z. If Vi = v this construction boils down to the one 



described in Section 1. 

It is straightforward to check that all constructions and results obtained in Sections 2- 
4 remain valid in this more general setting except only one point. In the definition of the 
proper pair of configurations one changes that the open segment between two mutually 
paired particles belonging to the interval [zi, Zi+i] for some i £ Z cannot exceed Vi (instead 
of v in the original definition). 

Therefore using the same arguments as in the proof of Theorem 1 we get its general- 
ization. 

Theorem 2 (Fundamental Diagram) Let z, Tj G {1,2, ...,r}, Vj G (0, v] be given and 
let p(x),p(z) be well defined and positive. Then 

V(x) = mm{l/p(z), l/p(x)}. (6.1) 



7 Discussion and generalizations 

1. The density of a configuration in the way how it was defined in Section 2 depends 
sensitively on the statistics of both left and right tails of the configuration. A close look 
shows that in fact if all particles move in the same direction, say right, one needs only 
the information about the corresponding (right) tail, which allows to expand significantly 
the set of configurations having densities and for which our results can be applied. 

For a configuration x G X by a one-sided particle density we mean the limit 

p(x) : = lim p(x, [0,£}) 
The upper an lower one-sided densities correspond to the upper and lower limits. 

Theorem 3 All previous results formulated in terms of "two-sided" densities remain valid 
if one replaces the usual particle density p to the one-sided density p. 

Proof. The key observation here is that the movement of a given particle in a config- 
uration x l G X depends only on particles with larger indices. Therefore if one changes 
positions of all particles with negative indices the particles with positive indices will still 
have the same average velocity. On the other hand, by Lemma 3 the average velocity 
does not depend on the particle index. This allows to apply the following trick. 

For each configuration x G X of density p(x) we associate a new configuration x G X 
defined by the relation: 

fx • if % > 

X' '.— \ 

\xq + i/p(x) otherwise . 

Then obviously p(x) = p(x) = p(x). 

Therefore for all purposes related to the average velocities all results valid for the 
configuration x remain valid for x as well. □ 

2. A close look to the proof of Theorem 1 may lead to the hypothesis that for a given 
configuration of obstacles z "typical" trajectories eventually become supported only by 
the locations belonging to the extended configuration z. Indeed, the calculation of the 
average velocity is given exactly for the configurations of this type. Let us show that this 



(7.1) 



is not the case even in the simplest setting when the configuration of obstacles has density 
one and is supported by integer points (i.e. a single obstacle at each point of Z). Let 
p(x) > 2 and let the initial configuration has at least 2 particles at each half-integer point 
\ Z. Then for each v > 1/2 and t > 1 the configuration x 1 is supported by the lattice 
ZU \ Z. 

3. Question about the existence of p(z). Let p(z) > be well defined. For each given v > 
we say that the configuration z is regular if the density of the corresponding extended 
configuration z is well defined. In the topology induced by the uniform metric in the space 
of sequences there exists an open set containing irregular configurations z. Nevertheless 
assuming a reasonable model for the creation of the configuration of obstacles one can 
show that a "typical" configuration of obstacles is regular. 

Indeed, assume that the sequence z is a realization of an ergodic stochastic process 
with stationary increments. Then Birkhoff Ergodic Theorem implies that Cesaro means 
for the sequence of fractional parts {(z i+ i — Zi)/v} converge to the limit and thus our 
claim follows. 

If p+{z) 7^ P-{z) then one may consider upper/lower particle velocities and the corre- 
sponding statement for the Fundamental Diagram may be rewritten as follows. 

Theorem 4 Let v,p(z) > be given. Then 

V±{x) =mm{l/p T {z), l/p T {x)}. 

The proof of this result follows basically the same arguments as in the case of Theo- 
rem 1 but some additional technical estimates related to partial limits are necessary. 

4. Assume now that we consider our original setting (i.e. Tj = Vj G Z) and the local 
velocity of the i-th particle does not depends on space (as in Section 6) but depends on 
time according to a given collection of local velocities {vj}, i.e. v\ stands for the local 
velocity of the i-th particle at time t. Thus using the notation introduced in Section 1 we 
get := min (A;(x, z), vfj and := x\ + £|. 

Theorem 5 For any given configuration of obstacles z such that Zk — > oo and two one- 
particle configurations x := {x } ^ x := {xo} there exists a sequence of local velocities 
{vq} such that \xq — x f \ > Ct for all t G Z + and some C > 0. Hence average particle 
velocities cannot coincide. 

The proof is based on the observation that one can choose {v^} such that each time 
t when the x-particle meets an obstacle it makes a "full" step equal to {v^} while the 
f-particle which meets with a different obstacle makes only a half step. 

This result shows that a direct generalization of our pure deterministic setting to the 
random case is not available. 



References 

[1] Angel O. The Stationary Measure of a 2-type Totally Asymmetric Exclusion Process, 
J. Combin. Theory Ser. A, 113:4(2006), 625-635. [math.CO/0501005] 

[2] Blank M. Ergodic properties of a simple deterministic traffic flow model. J. Stat. 
Phys., 111:3-4(2003), 903-930. [math/0206194] 

[3] Blank M. Hysteresis phenomenon in deterministic traffic flows. J. Stat. Phys. 120: 
3-4(2005), 627-658. [math.DS/0408240] 

[4] Blank M. Metric properties of discrete time exclusion type processes in continuum. 
J. Stat. Phys., 140:1(2010), 170-197. [math.DS+math.PR/0904.4585] 

[5] Borodin A., Ferrari P.L., Sasamoto T. Large time asymptotics of growth models on 
space-like paths II: PNG and parallel TASEP. [math-ph/0707.4207 ] 

[6] Comtet A., Majumdar S.N., Ouvry S. and Sabhapandit S. Integer partitions and 
exclusion statistics: limit shapes and the largest parts of Young diagrams. J. Stat. 
Mech. (2007) P10001. [cond-mat.stat-mech/0707.2312] 

[7] Evans M. R., Rajewsky N., Speer E. R. Exact solution of a cellular automaton for 
traffic. J. Stat. Phys. 95(1999), 45-98. 

[8] Evans M.R., Ferrari P.A., Mallick K. Matrix representation of the station- 
ary measure for the multispecies TASEP. J. Stat. Phys., 135:2(2009), 217-239. 
[math.PR/0807.0327] 

[9] Gray L., Griffeath D. The ergodic theory of traffic jams. J. Stat. Phys., 105:3/4 
(2001), 413-452. 

[10] Ispolatov I., Krapivsky P.L. Particle systems with stochastic passing. Phys. Rev. E 
61(2000), R2163-R2167. [cond-mat/9912293] 

[11] Kornfeld E.P., Sinai Ya.G., Fomin S.V. Ergodic theory. M.: Nauka, 1980. (Springer 
Verlag, 1982) 

[12] Krug J., Garcia J. Asymmetric particle systems on R. J. Stat. Phys., 99:1-2(2000), 
31-55. [cond-mat.stat-mech/9909034] 

[13] Liggett T.M. Interacting particle systems. Springer- Verlag, NY, 1985. 

[14] Nagel K., Schreckenberg M. A cellular automaton model for freeway traffic, J. 
Physique I, 2 (1992), 2221-2229. 

[15] Penrose M.D. Existence and spatial limit theorems for lattice and continuum particle 
systems. Probab. Surveys, 5 (2008), 1-36. 

[16] Zielen F., Schadschneider A. Exact Mean-Field Solutions of the Asymmetric Ran- 
dom Average Process. J. Stat. Phys., 106:1-2 (2002), 173-185 [cond-mat.stat- 
mech/0 104298] 



